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CERN collaboration NA50 has measured charmonium and Drell-Yan dimuon production in Pb+Pb 
collisions. Parton scattering broadens the transverse momentum, p T , distributions for these pro- 
cesses. We predict that (p^,) will flatten in Pb+Pb collisions as a function of the neutral transverse 
energy of hadrons, Et, in contrast to the almost-linear rise seen in S+U — > ip + X. If seen, such a 
flattening will support hadronic explanations of charmonium suppression. 
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The NA50 collaboration at the SPS has reported a sup- 
pression of ip production in Pb+Pb collisions relative to 
Drell-Yan dimuon production as the neutral transverse 
energy of hadrons, Et, is increased Jjj]. They further pre- 
sented a striking 'threshold effect' by comparing the data 
to S+U results as a function of a calculated quantity, the 
mean path length of the ip through nuclear matter, L. 
The NA50 comparison is shown in fig. 1; the open sym- 
bols are from ref. pi. In this note, we demonstrate that 
the relation between L and the measured Et is model 
dependent. However, we point out that measurements of 
the centrality dependence of the transverse momenta of 
Drell-Yan dimuons essentially provide an experimental 
determination of L. We argue that such a determination 
can provide vital evidence for - or against - the threshold 
behavior, which has been linked to the onset of quark- 
gluon plasma formation in the Pb system Jl|,||. 

To see how the path length affects the 'threshold' inter- 
pretation of the NA50 comparison, we replot the Pb+Pb 
data using L(Et) calculated with eq. (||) and the realistic 
nuclear densities (|] that were employed in Q. We see 
that the appearance of the data is very sensitive to small 
changes in the definition of L. With the realistic L, one 
no longer gets the impression that the Pb+Pb data "de- 
parts from a universal curve." This departure has been 
cited as evidence Q that Pb collisions cross a threshold 
and form quark-gluon plasma. 

In ref. Q, we observed that the Pb+Pb data are in 
good accord with predictions of charmonium suppression 
using a hadronic comover model m. While these pre- 
dictions are supported by cascade calculations [Q, the 
plasma explanation cannot be excluded. We attributed 
the behavior of fig. 1 to a geometric saturation that oc- 
curs in the symmetric Pb+Pb system, but not in asym- 
metric S+U collisions; we discuss this below. We see in 
fig. 1 that the saturation phenomena, though softened, 
does not vanish with our calculated L(Et)- Saturation - 
the flattening of L{Et) - indeed occurs in our scenario, 
but only as a consequence of geometry. 

Important information on the suppression mechanism 
can be extracted from the nuclear dependence of the ^'s 



transverse momentum, p T . The hadronic suppression in 
ref. [QJ^] does not modify the transverse momentum de- 
pendence of ip production appreciably ||. To account 
for the nuclear modification measured in pA and S+U 
collisions, initial-state parton scattering was introduced 
l^tTof. Bodwin, Brodsky and Lepage and, independently, 
Michael and Wilk had predicted that initial- and final- 
state parton scattering can modify the momentum dis- 
tributions of hard interactions in nuclear collisions ]Tl| . 
Evidence for such scattering has been seen in a variety 
of hadron-nucleus experiments I2-I4[. 
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FIG. 1. Comparison of NA50 Pb+Pb (open circles) and 
NA38 S+U (open triangles) data as in W. Data replotted 
with a realistic L(Et) from eqs. (^[0) appear different. The 
gray area represents the uncertainty in this presentation of 
the data due to L(Et)- 

Initial-state scattering JIlJ is elastic at the parton 
level and broadens the p T distributions of charmonium 
and Drell-Yan production without affecting the p T - 
integrated yields ||[l0|]. Partons essentially undergo a 
random walk in momentum space, so that (p~,) grows lin- 
early with L. Initial-state scattering provides the only 
nuclear effect proposed so far that can modify the p T dis- 
tribution in the Drell-Yan process in a manner consistent 
with data |l2|,|f3|. While other effects such as comover 
scattering can alter the p T distribution in charmonium 
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production, we find in practice that the contributions to 
(p 2 ) are small if the comovers are assumed to be hadrons 
(see || and the discussion below). Measurements of (p 2 ) 
- particularly for dimuons - can therefore provide exper- 
imental information on L. 

If (p 2 ) increases with L as in hadronic models of char- 
monium suppression incorporating initial-state scatter- 
ing, then the geometric saturation in Pb collisions im- 
plies a flattening of (p 2 ). In contrast, if the thresh- 
old interpretation JlJ is correct and the Pb results are 
the consequence of a new contribution to ip suppression 
from quark-gluon plasma, then one would not expect 
(p 2 ) to flatten with Et- If observed, a flattening would 
lend strong support to the comover explanation. Alter- 
natively, if no flattening is observed, then our hadronic 
model would be excluded. Plasma suppression is cur- 
rently thought to be strongly p T dependent jl5| , although 
models with a weaker dependence (and without thresh- 
olds) have been presented Jl(J. 

Plasma issues aside, measurements of(p 2 ) can help re- 
move the model dependence of fig. 1 |Lj. Gerschel and 
Hiifner (l7j pointed out that a plot similar to fig. 1 but 
containing only measured quantities can be constructed 
by replacing L with the measured (p 2 ) of the ip. Mea- 
surements of the (p 2 } of the Drell-Yan process, though 
more difficult, can provide a clearer determination of L. 

In this paper we update the description of initial-state 
scattering from ref. |^] using more recent information 
from FNAL pA experiment E772 |D| . We then derive an 
expression for the path length L, and study its behavior 
as a function of Et for S+U and Pb+Pb collisions. Fi- 
nally, we combine these results to predict the behavior of 
(p 2 ) for Pb+Pb collisions. 

In a hadron-nucleus collision, a parton from the pro- 
jectile can suffer soft quasielastic interactions as it crosses 
the nuclear target on its way to the hard process. We fol- 
low |9| and assume that in each nucleon-nucleon, NN, 
subcollision there is a fixed probability </> that the parton 
is affected. In a hadron-nucleus collision, the (p 2 ) of the 
dimuon or ip is then increased by: 



Ap 2 ^(p 2 )-(p 2 ), 



X 2 (n A - 1), 



(1) 



where n A is the number of NN subcollisions that the 
projectile suffers in the target and A 2 oc <p determines the 
increment to the projectile parton's p 2 from each subcol- 
lision. We subtract '1' to eliminate the hard subcollision 
that produces the ip or dimuon from n A and introduce 
(Pt) nn to account for the ^-independent contribution 
from that subcollision. The impact parameter averaged 
total number of subcollisions grows with the target ra- 
dius R A ps 1.2 A 1 / 3 as n A » 3a NN p Q R A /2 rj 0.77 A 1 / 3 , 
where po is the average nuclear density and a NN ~ 32 mb 
is the inelastic NN cross section. 

In ref. transverse momentum distributions mea- 
sured in hadron-nucleus collisions by CERN NA10 and 
NA3 [pi were used to fix {p 2 T ) NN and A for Drell-Yan 
and ip production. In fig. 2, we compare Ap 2 from 



cq. (|l|) with more recent data from a larger number of nu- 
clear targets presented by FNAL E772 Jl3|] . The random 
walk behavior, Ap 2 oc A 1 / 3 , is evident. This approx- 
imation agrees roughly with calculations using realistic 
nuclear densities as in ref. Jig] . In this work we take 
A M+M - pb 0.18 ± 0.01 GeV and A^ « 0.36 ± 0.03 GeV for 
Drell-Yan and ip production. These values are chosen to 
agree with E772 and NA3 data, respectively. They are 
compatible with the idea that initial-state scattering is 
a soft process. We will use NA38 S+U data to obtain 
(P 2 T ) nn below. 

Observe that our Drell-Yan value is somewhat smaller 
than the 0.24 ± 0.05 GeV value used in ref. g. This 
spread perhaps reflects the systematic uncertainty (not 
included in the fit uncertainty) in comparing experiments 
with different kinematic coverages. E772 did not report 
Ap 2 for ip, but did give results for T production, a pro- 
cess that likely has similar initial-state interactions. A 
fit to that data implies Ax ~ 0.46 GeV, which is larger 
than our A^ = 0.36 GeV. The values that we have taken 
above are the smaller of the choices. 
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FIG. 2. Calculated Ap 2 T from eq. ^ for n A « 0.77 A 1/z 
compared to FNAL E772 data on the Drell-Yan process and 
T production Jl3|. Diamonds are calculated using (^]) and 
JVV nuclear densities from 

We remark on the interpretation of the parameters in 
eq. (|I]). In ref. ||, the value of the ratio (A^,/A /J + /i - ) 2 ~ 
2.3 had been attributed to the 9/4 enhancement of the 
cross section for gluon-gluon relative to quark-gluon for- 
ward scattering. The initial-state partons are predom- 
inantly gluons for -0 production and quarks and anti- 
quarks for the Drell-Yan process, because the primary 
production mechanisms for ip and dilepton production 
are gluon fusion and quark-antiquark annihilation, re- 
spectively. In contrast, the parameters that we now use 
imply a ratio (A^/A^^- ) 2 ~ 3.9 that is ~ 70% larger 
than the earlier estimate, while E772 values alone sug- 



gest (At/a„+„- y 



6.8. The perturbative estimate 



2 



(A^,/A M + M -) 2 ~ 9/4 |lC|]li| may not be correct in detail, 
since A 2 involves small momentum transfers. Neverthe- 
less, this disagreement is quite large. 

If taken literally, a larger ratio can suggest that the ip 
also undergoes final-state elastic scattering as it escapes 
the nucleus. During its escape, the nascent ip is not a 
fully formed hadron but, rather, a cc pair. This pair can 
be in a color octet state and, therefore, can scatter essen- 
tially as a gluon does. Fortunately, this final-state octet 
scattering also follows eq. ([[]), so that we can describe 
this effect in conjunction with initial-state gluon scatter- 
ing by replacing A 2 with 2A 2 . Together with the 9/4 color 
factor, this factor can roughly account for the ratio of A 2 
in ip and /i + /J~. However, this observation is speculative. 
Therefore, we will not include the factor of two explicitly, 
but merely stress the possibility that the empirical A for 
charmonium production receives contributions from both 
initial- and final-state octet scattering. 
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FIG. 3. NA50 L(E T ) [1] (points) compared to calcula- 
tions for the JVV nuclear densities Q (solid) and for a 
sharp-surface approximation (dot-dashed). 

In a nucleus-nucleus collision, both projectile and tar- 
get partons scatter. We then write 

Apl^(pl)-{pl) NN =X 2 (n A + n B -2). (2) 

The NA50 path length is 

L = (n A + n B )/2cr NN p Q = n/2a NN p . (3) 

The relation between n and the impact parameter b de- 
pends on the collision geometry. A nucleon at a trans- 
verse position s in the projectile can undergo 



n A 



dz' p A (s,z'), 



subcollisions prior to its hard interaction 
nal position z, where p A is the projectile 
expression for the number of subcollisions 
nucleon at a transverse position b — s in 
similar. Since the hard process occurs with 
density cx p A (s, z)ps(b — s,z'), the averaj 
subcollisions for a given b is 



(4) 

at longitudi- 
density. The 
suffered by a 
the target is 
a probability 
ie number of 



n(b) 



Tai 



d 2 sT A (s)T B (\b-s\)[T A (s) 



+T B (\b-s\)}, (5) 

where T Ai b — X^^ 2 ' Pa, b are the nuclear thickness 

functions and T A b = [d 2 s T A (s)TB{\b — s|). An alterna- 
tive derivation of eq. (|) in the spirit of rcf. is given 
in ref. @. 
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FIG. 4. Average impact parameter as a function of Et for 
realistic nuclear densities (solid) and for a sharp-surface ap- 
proximation (dot-dashed). NA50 L(Et) [1] (points) are in- 
cluded for comparison. 

To obtain the Et dependence of n - and therefore 
L(E T ) - we fold eq. (||) with the probability P(E T ,b) 
that a collision at impact parameter b produces trans- 
verse energy Et- This probability is related to the 
minimum-bias distribution by 



l (E T ) = Jd 2 bP{E T ,b). 



(0) 



The distribution P(ET,b) depends both on the collision 
geometry and the calorimeter type. We use the phe- 
nomenological P(Ex,b) relevant to the NA50 Pb+Pb 
collisions from ref. 0] ; this parametrization schematically 
incorporates fluctuations at fixed b due to both collision 
dynamics and detector effects. The average number of 
collisions for a given Et is then 



n(Ej 



,(T A (s)+T B (\b-s\)), 



(7) 



where we now average over the weighted product of the 
densities, 

(...) fd 2 bd 2 S P(E T ,b)T A ( S )T B (\b-s\)(...), 



(8) 

with N = Jd 2 bP(E T ,b)T AB (c.f. eq. (§V). 

We now compute L(Et) using eqs. ^,@) and (|J). In 
fig. 3 we compare the NA50 L{Et) Q to the path length 
calculated using two assumptions for the nuclear den- 
sity profile. Following ref. Q, we use the realistic three- 
parameter Woods-Saxon densities from deJager, deVries 
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and deVries (JVV) ||. We compared this to a sharp- 
surface approximation p = p Q(Ra — r). NA38 Q ob- 
tained L for S+U from ip p T data using a phenomeno- 
logical procedure that is essentially equivalent to eqs. (||, 
||), while NA50 calculated L assuming the sharp-surface 
approximation . Consequently, we see that the NA50 
Pb+Pb values agree with our sharp-surface results, while 
the NA38 S+U values are nearer to the realistic-density 
computations. [The calculation of L is not described in 
ref. so the agreement of our sharp-surface L(Et) 
with the NA50 Pb result is an important check.] Note 
that to compare to the NA50 calculations, we plot PqL 
rather than L because the authors of ref. [Q use a value 
pa = 0.138 fm- 3 . 

We see that L(Et) saturates in Pb+Pb collisions, as 
observed in ref. W . To see why saturation occurs in this 
system but not in S+U, we compare the NA50 L(Et) @] 
to the average impact parameter (b)(Er) in fig. 4. We 
use the JVV densities to compute (b) — (oTab) I (Tab) ■ 
[Note that NA50 reports similar values of (b)(E T ) 
For all but the highest E T bin addressed by the S+U 
measurement, we see that (b) is near ~ i?g = 3.6 fm or 
larger. In this range, increasing b dramatically reduces 
the collision volume and, consequently, L. In contrast, in 
Pb+Pb collisions (b) < R Ph = 6.6 fm for E T > 50 GeV, 
so that L does not vary appreciably. 
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FIG. 5. The calculated Et dependence of (p^,) for ip (up- 
per plots) and Drell-Yan dimuons (lower plots) in S+U and 
Pb+Pb collisions. The calculations are performed for real- 
istic nuclear densities (solid) and for the sharp-surface ap- 
proximation (dot-dashed). NA38 S+U data is from [ B2fl . 
The thick curves are computed for our preferred values, 



Av, = 0.36 GeV and X 



— 0.18 GeV, while the thin curves 



are for X^p = 0.46 GeV and X 



= 0.24 GeV. 



We now use eqs. S 0) to compute Ap 2 for both the 
Drell-Yan process and ip production as functions of Et 
in S+U and Pb+Pb collisions. Our results are shown 



in fig. 5. Parton-scattering calculations agree with S+U 
data from rcf. p^] . Observe that to compare to the NA38 
S+U data, we fit the data to extract (p 2 ) NN ~ 0.92 and 
1.07 GeV 2 for Drell-Yan pairs and ip respectively. Note 
that the ip value is somewhat smaller than the 1.23 ± 
0.05 GeV 2 reported by NA3 and used in ref. M; this is 
likely another indication of systematic uncertainties. 

To illustrate the spread in the Pb predictions implied 
by the uncertainty in the parameters, we compute the 



thin curves in Fig. 5 by taking A 



Ai 



from E772 and A^+p 



0.46 GeV 

0.24 GeV from |9j. With these 
values, we obtain the best agreement with the NA38 data 
for (p 2 T ) NN w 0.81 and 0.6 GeV 2 for dimuons and ip, re- 
spectively. The modified ip results are inconsistent with 
the S+U data. The modified dimuon results are consis- 
tent, although the agreement is somewhat better for our 
preferred value, A^+ M - = 0.18 GeV. 

The thick curves for Pb+Pb collisions in fig. 5 rep- 
resent our predictions. We expect (p 2 ) to increase by 
12.3% for ip as Et increases from 50 to 150 GeV. This 
represents a flattening of (p 2 in comparison to the 

S+U — > ip+X data, which show an 18.5% increase as Et 
varies from 30 to 90 GeV. The flattening would be even 
more dramatic if the NA50 sharp-surface approximation 
(the dot-dashed curve in fig. 5) were true — (p 2 ) would 
increase by only 6.8% for ip and 2.7% for dimuons. While 
such small differences seem difficult to resolve, agreement 
with S+U — » ip + X data is better for the realistic L 
than for the sharp-surface result. Furthermore, NA38 
po| , pT| used this ip data to obtain L values in agreement 
with our realistic calculation and quite distinct from our 
sharp-surface result. However, we stress that the most 
direct extraction of L comes not from ip but from Drell- 
Yan dimuons. To check that the NA50 plot Q is correct, 
i.e. to decide between the open circles and filled circles in 
fig. 1, Drell-Yan data in Pb+Pb would have to establish 
a 2.7% increase in (p 2 ). This requires dimuon data far 
more precise than that for S+U collisions. 

One can ask if elastic scattering of the ip by co- 
movers can affect its transverse momentum distribution. 
As mentioned earlier, comover absorption has only a 
marginal effect on the p T distribution; its influence on 
(p 2 T ) is negligible However, elastic comover scatter- 

ing can in principle allow the ip to be pushed transversely, 
adding to its p T . The argument [p|Jl5| that elastic scat- 
tering with hadronic comovers is negligible follows from 
photoproduction data. The measured elastic ipN cross 
section is ~ 0.079 + 0.012 mb p3| , accounting for only 
~ 2 — 4% of the total cross section. The corresponding 
mean-free path for elastic scattering with hadronic co- 
movers greatly exceeds the estimated size of the system, 
so that the ip will not follow the comover flow. The gen- 
eral relation between scattering and flow is discussed in 
ref. [g4|. One possible hole in this argument is that ipir 
scattering in the comover gas is likely below the DD dis- 
sociation threshold and, therefore, predominantly elas- 
tic. Another possible hole is that comovers need not be 
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hadrons. 

The consistency of the NA38 S+U -> V + X data in 
fig. 5 with (0) supports our neglect of elastic comover 
scattering. Nevertheless, if Pb+Pb experiments find a 
(p 2 ) for tp that is larger than our prediction and if other 
hadronic species show evidence of substantial transverse 
flow, it will be necessary to introduce elastic scattering 
into the comover scenario. This is best done in the con- 
text of cascade models. 

In summary, we have predicted the nuclear enhance- 
ment of (p 2 T ) in Pb+Pb collisions as a function of Et for 
charmonia and Drell-Yan dimuons, assuming that this 
enhancement is caused by quasielastic parton scattering. 
Such scattering has been included P, ^0|Jl5| in hadronic 
models of tp suppression, where it is essential for describ- 
ing the p T dependence of pA, O+U and S+U data. We 
stress that the Drell-Yan process is unaffected by final- 
state comover or plasma interactions, so that a compar- 
ison of both ip and calculations to data can disen- 
tangle model uncertainties in L(Et) from new physics. In 
particular, the L extracted from dimuon measurements 
can confirm or disprove the threshold behavior claimed 
by NA50. 

The parameters in the model were revised from an ear- 
lier work H to describe the latest E772 pA — > fi + fi" + X 
data. Our revision implies a ratio (A^,/A Al + A1 - ) 2 ~ 3.9, 
larger than that extracted earlier || , perhaps indicating 
that additional final-state scattering of the octet cc oc- 
curs in the charmonium case. More precise pA — > ip + X 
measurements are needed to explore this very interesting 
possibility. 
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